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ABSTRACT 



IE 1207.4— 5209, the peculiar Central Compact object in the G296. 5+10.0 supernova remnant, has been proposed to 
be an "anti-magnetar" - a young neutron star born with a weak dipole field. Accretion, possibly of supernova fallback 
material, has also been invoked to explain a large surface temperature anisotropy as well as the generation of peculiar 
cyclotron absorption features superimposed to its thermal spectrum. Interestingly enough, a faint optical/infrared 
source was proposed as a possible counterpart to IE 1207.4—5209, but later questioned, based on coarse positional 
coincidence. Considering the large offset of IE 1207.4— 5209 with respect to the center of its host supernova remnant, 
the source should move at ~ 70 mas yr _1 . Thus, we tested the association by measuring the proper motion of the 
proposed optical counterpart. Using Hubble Space Telescope (HST) observations spanning 3.75 years, we computed a 
3cr upper limit of 7 mas yr _1 . Absolute astrometry on the same HST data set also places the optical source significantly 
off the 99% confidence Chandra position. This allows us to safely rule out the association. Using the HST data set, 
coupled to ground-based observations collected at the ESO/ Very Large Telescope (VLT), we set the deepest limits ever 
obtained to the optical/infrared emission from IE 1207.4—5209. By combining such limits to the constraints derived 
from X-ray timing, we rule out accretion as the source of the thermal anisotropy of the neutron star. 

Key words. Stars: neutron, Stars: individual: IE 1207.4—5209 



1. Introduction 

The X-ray source IE 1207.4—5209 was discovered with the 
Einstein satellite (Hclfand & Becker 1984) close to the 
center of G296. 5+10.0, a ~ 7 kyr old s upernova remnant 
(SNR) located at a distance of ~ 2 kpc (Roger et al.||1988 



peculiar binary system ( ZavlinetaL 12004 



2007). However, more recently, Gotthelf & Halpernl (12007) 
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Woods et al. 



Giacani et al.pOOO). It was the second thermally-emitting 



radio-quiet, Isolated Neutron Star (INS) candid ate found 
inside a SNR, after IE 161348-5055 in RCW 103 flTuohyfe 
Garmire|1980 |. A handful of similar sources, discovered in- 
side young SNRs , are dubbed, as a class, Central Co mpact 
Objects (CCOs, |Pavlov et al.||2002| |de Luca||2008[ for a 



recent review) . Although their properties are not well un- 
derstood, CCOs are supposed to be the youngest members 
of the radio-quiet INSs family. 

IE 1207.4—5209 is one of the most peculiar and most 
observed Galactic X-ray sources. Pulsation s at 424 ms 
were discovered with the Chandra satellite ( Zavlin et al. 



provided robust evidence that IE 1207.4—5209 is a very 
stable rotator. The upper limit to its period derivative 
(P < 2.5 x 10~ 16 s s _1 at 2a) yields an INS carachteristic 
age r c > 27 Myr, exceeding by 3 orders of magnitude the 
age of the SNR, a rotational energy loss E < 1.3 x 10 32 erg 
s _1 , and a very small dipole magnetic field, B< 3.3 x 10 11 
G. Thus, IE 1207.4—5209 could be a young, weakly magne- 
tized INS, born with a spin period very similar to the cur- 
rent one. Evidence for similar, low magnetic fields has been 
obtained for two other members of the CCO class, namely 
CXOU J185238 6+00 4020 at the center of the Kes 79 SNR 
ffialpern et"alj|2007[) and RX J0822-4300 in Puppis A 
(Gotthelf & Halpcrn 2009), adding support to the scenario 
of CCOs as "anti-magnetars" (with the remarkable excep- 
tion of the puzzling source in RCW 103, De Luca et al 



2000), proving the source to be an INS. Early timing in- 
vestigations hinted at a non-monotonic period evolution 
of IE 1207.4—5209 suggesting that the source could be a 

Send offprint requests to: deluca@lambrate.inaf.it 
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2006) 



What makes IE 1207.4-5209 unique among all INSs is 
its X-ray spectrum. Three (possibly four) broad absorption 
features - at regularly spaced energies (0.7, 1.4, 2.1 and pos- 
sibly 2.8 ke V) - are visible over a the r mal two-temperature 



continuum (Mereghett i et al. 2002| jSanw al et al. 
iDe Luca et al 



[2002 
The rJeptrTo 



|Bignami et al . 2003; L>e Luca etaT N2004J 

such features varies as a function of the rot ational phase 

(|Mereghetti et al.||2002) |De Luca et al.||2004[). The nature 



A. De Luca et al.: HST and VLT observations of IE 1207.4-5209 



of the spectral features of IE 1207.4—5209 has been debated 
since their discovery, possible interpretations being atomic 
transition lines in the NS atmosphere or cyclotron features 



in the plasma surroun ding the star (Sanwal et al. 2002 



Mereghetti et al.|2 002). The latter interpretation, favourec 



by the harmonic energy spacing of the features (Bignami 



et al. 2003 De Luca et al. 



the emerging picture of IE 



20041), is fully consistent with 
207.4-5209 as a weakly mag- 



netized neutron star: assuming the 0.7 keV feature to be 
the fundamental electron cyclotron line yields a measure of 
the magnetic field of 8 x 10 10 G, i.e. below the value derived 
from the upper limit on the pulsar spin down. 

Many puz zles remain to be solved for IE 1207.4—5209. 
As stated by |Gotthelf k Halpern| ( |2007[ ), it is difficult to 
explain the luminosity and temperature of the observed hot 
thermal spectral component within the frame of a weakly 
magnetized INS. Moreover, the physics underlying the pe- 
culiar absorption features (most likely due to cyclotron pro- 
cesses, Liu et al. 2006 |Suleimanov et aL"1|2010| |Potekhin 



2010) has to be understood. The simplified model proposed 



By" Liu et al. (2006) requires a high (and steady) electron 
density in the NS magnetosphere above the polar cap. Low- 
level accretion, possibly of supernova fallback material, has 
been invoked to ease the problem in both cases. This would 
point to the existence of a debris disk surrounding the INS 
(which could be detected in the optical- infrared range), a 
long sought astrophysical object, so far possibly observed 
only in th e case of the Anom alous X-ray Pulsar (AXP) 4U 
0142+61 flWanget al.|[2006| ). 

Very deep imaging of the field of IE 1207.4-5209 have 
been performed both from the ground with the ESO 
Very Large Telescope (VLT) and with the Hubble Space 
Telescope (HST). Optical VLT observations (De Luca et 
al. 2004) did not reveal any potential counterpart down 
to R~27.1 and V^27.3 while observations in the optical 
with the HST and in the near infrared (NIR) with the VLT 
showed the presence of a faint source (hereafter "source 
Z" ) close to the Chandra X-ray position, with magnitudes 



m F555W ~ 26.4 and Ks~20.7 (Pavlov et al. 



et al.|2006| ). However, the ass ociation" with IE WM - 5209 



2004; 



Fcscn 



was soon after questioned by Mignani et al. ( 2007a ) on the 
basis of precise absolute astrometry of the HST images, 
which showed a positional offset of source Z with respect 
to the Chandra coordinates. The same source was observed 



in the NIR by |Wang et all ( 2007 1 , who reported very red 

M dwarf, and also questioned its 



colours, consistent with an 

possible association to IE 1207.4 

consistency with the Chandra position. Wang et al. (2007) 



-5209 because of the in- 



also observed the field with Spitzer at 4.5/mi and at 8.0/mi 
but did not detect any source at the target position. 

Here we report on a different, independent test to assess 
the association of source Z to IE 1207.4—5209, using multi- 
epoch data collected with the HST (Sect. ^. The same HST 
dataset, completed by ground based data collected with 
the VLT, is also used to derive stringent constraints on the 
optical/infrared emission from IE 1207.4-5209 (Sect. [3). 
Results are discussed in Sect. |4j 



2. Association of source Z to IE 1207.4-5209: an 
HST test. 

G296. 5 +10.0 has a remarkable, well defined b ilateral sym- 
metry (|Roger et ai]|1988| |Storey et ai]|1992[). Very likely, 



the explosion site lies on the symmetry axis of the SNR, 
but the current position of IE 1207.4—5209 is significantly 
offset from the apparent center of the host SNR. Indeed, 
the ge ometrical center position evaluated by |Roger et al. 
(1988 ) is ~ 8' to the south west of the X-ray source (De 
Luca et al~||2004 1 . Assuming for the system an age of 7 000 
years, such a displacement would imply a proper motion 
of ~ 70 mas yr _1 , corresponding to a projected velocity of 
~ 640 km s _1 , consistent with the observed velocity dis- 
tribution for radio pulsars (Hobbs et al. 2005). This offers 



a natural way to test the putative identification: if indeed 
associated with IE 1207.4—5209, source Z should show a 
significant proper motion. Thus, we used multi-epoch HST 
observations to search for an angular displacement of source 
Z. 



2.1. HST observations 

We observed the field of IE 1207.4-5209 with the HST 
on May 8th 2007 (Programme 10791). Our observa- 
tions were originally scheduled for execution with the 
Wide Field Channel ( WFC) of the Advanced Camera for 



Surveys (ACS) ( |Clampin et al. |2000| [Sirianni et al.||2005[ ) 



Unfortunately, the ACS/ WFC was put in idle state on 
January 2007 due to a failure of the on-board electron- 
ics. Our observations were re-scheduled and executed with 
the Wide Field Planetary Camera 2 (WFPC2). A set of 
four 500s exposures were obtained during one spacecraft 
orbit, through the 814W filter (A = 8012A; AA = 1539 
A). In order to exploit the maximum spatial resolution 
for the proper motion measurement, IE 1207.4—5209 was 
placed at the centre of the Planetary Camera (PC) chip 
(O'.'045/pixel). 

Our new data add up to observations collected with 
the ACS' on July 28th and August 7th 2003 (Programme 
9872) and available in the public HST archives . This first- 
epoch HST dataset allowed Pavlov et al. ( 2004 ) to pick up 



Source Z as a possible counterpart to IE 1207.4—5209. The 
WFC (O'.'050/pixel) was used in both visits. Two sequences 
of 4 and 5 exposures were obtained through the broad- 
band filters 555W (A = 5346A; A A = 1193 A) and 814W 
(A = 8333A; A A = 2511 A), for a total integration time 
of 12800 s and 10200 s, respectively. The complete dataset 
spans a time baseline of ~ 3.75 years. 

We downloaded the data from the Space Telescope 
European Coordinating Facility (ST-ECF) Science Data 
ArchiveR On-the-fly data reduction (bias and flat-field 
correction) and flux calibration were applied using the 
Space Telescope Science Data Analysis Software (STSDAS) 
through the ST-ECF Data Archive pipeline. To filter cosmic 
ray hits, single WFPC2 exposures were combined and av- 
eraged using the STSDAS task combine, while single ACS 
exposures were combined using multidrizzle which also 
produces a mosaic image of the two ACS chips and applies 
the correction for the geometric distortions of the camera. 

2.2. Relative astrometry 

We have used the algorithm that we successfully applied 
in a series of previous proper motion measurements with 
HST (see e.g. |De Luca et al.|[2000| |Caraveo et al.||2001" 



Mignani et al. |2002[ |De Luca et all, 2007[ |Mignani et al 



.stecf.org/archivo 
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Fig. 1. 16" x 11" cutout of the 2003 ACS image of 
IE 1207.4-5209 obtained through the 814W filter. The 
Chandra position of the target, as computed using the ACIS 
and the HRC observations, is shown. The circles (C'62 ra- 
dius) corresponds to the uncertainty of the CCO position 
computed from the error of the Chandra coordinates and 
the accuracy of the 2MASS absolute astrometry. North to 
the top, East to the left. 



20101. We used only the ACS F814W bandpass image in 
order to avoid introducing systematic effects in the com- 
parison of images taken through different filters. Briefly, a 
relative reference frame is defined by a grid of good (not ex- 
tended, with good signal to noise but not saturated, not too 
close to the CCD edges) reference sources common to all 
images and spread homogeneously across the common field 
of view. We selected 32 such sources and we computed their 
positions by fitting a 2-D gaussian profile to their brightness 
profile. The centering accuracy was between 0.02 and 0.07 
pixels per coordinate for the WFC and between 0.04 and 
0.08 pixels for the PC depending the source brightness. The 
position of source Z was evaluated in the same way, with 
an uncertainty of 0.08 and 0.04 pixels in the PC and in the 
WFC image respectively. Next, for the PC image, we cor- 
rected t he pixel coordinates for the effects of the "34th row 
effect" ( Anderson fc King|1999[) as well as for the geometric 
distortion ( |Anderson fc Kin g 2003). Then, we assumed the 
2007 image as a reference and we aligned the coordinate 
grid along Right Ascension and Declination using the well- 
measured HST roll angle (^ 0.1° uncertainty). We then 
registered the 2003 reference frame to the 2007 one by fit- 
ting a 5-parameter transformation (including x and y shifts 
and scale factors, as well as a rotation angle). We rejected 5 
objects yielding > 3<7 residuals using an iterative cr-clipping 
algorithm. We ended up with a r.m.s. residual of 0.17 PC 
pixels on the frame registration using 27 reference sources. 
Finally, we applied the transformation to the coordinates 
of source Z in order to evaluate its possible displacement 
between the two epochs. However, such a displacement was 
found to be much smaller than the overall uncertainty of 
0.19 PC pixels. Taking into account the PC plate scale 
(O'.'04554/pixel) as well as the time baseline covered by our 
dataset, we can set a 3cr upper limit of 7 mas yr _1 to the 
overall proper motion of source Z. This limit is much lower 
than the expected proper motion of ~ 70 mas yr _1 implied 
by the association between IE 1207.4-5209 and the SNR. 



We note that the proper motion test relies on assump- 
tions for the supernova remnant age and expansion center 
position. The age is quite uncertain. However, even assum- 
inga value 3 times higher than estimated by |Roger et al.| 
( 1988 ) , we would expect a proper motion of ~ 25 mas yr _1 . 



The good bilateral symmetry of G296. 5+10.0 suggests 
that the true expansion center should lie not far off the 
apparent symmetry axis. However, expansion of the super- 



nova remnan t could have been anisotropic. Indeed, |Roger 
|et al. (1988) report a smaller radius of curvature for the 
Eastern side of the SNR with respect to the Western side, 
suggesting an expansion center to the East of the geometri- 
cal center (i.e. closer to the CCO than the geometrical cen- 
ter), although they did not quantify such effect. We might 
assume that the true expansion center is placed halfway 
between IE 1207.4-5209 and the SNR geometrical center, 
yielding a factor 2 error in our estimate. Even in this case 
(and even combining the two pessimistic assumptions), the 
expected proper motion would have been larger than our 
upper limit. Thus, this strongly argues against association 
of source Z to IE 1207.4-5209. 



2.3. Absolute astrometry 

A second, model-independent, test on the possible associ- 
ation between IE 1207.4—5209 and source Z is obviously 
represented by the comparison of their positions. Before 
registering the IE 1207.4—5209 position on our HST images 
we decided to recompute anew its X-ray coordinates from 
the Chandra data, to account for possible corrections in the 
satellite aspect solution. To this aim, we retrieved all the 
archival Chandra imaging observations of IE 1207.4—5209. 
The first observation was collected with the Advanced CCD 
imaging Spectrometer (ACIS-S) on 2003, June 6 (Obs.Id. 
3913), totalling ~ 20 ks; the second one was collected with 
the High Resolution Camera (HRC-I) on 2003, December 
28 (Obs.Id. 459 3) and lasted ~ 50 k s. Data analysis was 
performed as in De Luca et al. (2009). We note that while 
Mignani et al. (2007) used only the coordinates derived 
from the ACIS-S data set, Wang et al. (2007) used co- 
ordinates derived from both data sets. The position of 
IE 1207.4-5209 turned out to be a( J2000)=12 /l 10 m 0.92 s , 



<5(J2000) = -52° 26'28"40 in the ACIS observation and 
a(J2000)=12' 1 10 m 0.88 s , 6(J2000)= -52° 26'28'.'61 in the 
HRC observations. The two positions agree within ff.'A, as 
expected according to Chandra astrometric accuracy (0'.'6 
at 90% confi dence levep and are consistent with those 
computed by |Wang et aT| ( |2007 ) using the same data. To 
evaluate the accuracy of the Chandra absolute astrome- 
try, we cross-correlated the position of X-ray sources de- 
tected within 3' of the optical axis with that of stars in 
the T wo Micron All Sky Survey (2MASS, |Skrutskie et al. 
2006 ) catalog. We found two matches between 2MASS stars 
and sources detected both by the ACIS and HRC. Based 
on the source density, we expect a chance alignment of 
one 2MASS source to a Chandra source within the nom- 
inal 90% error region to have a probability of 0.6% and 
of 0.2% for ACIS and HRC, respectively. Thus, it is very 
likely that the two 2MASS sources are the IR counterparts 
of the matched Chandra sources. The difference between 
the Chandra and the 2MASS coordinates of these sources 
is of ~ 0"4 and is consistent with the expected astromet- 



http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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ric accuracy of Chandra. Offsets along Right Ascension and 
Declination range from O'.'OT to 0'.'33 and have different di- 
rections. Thus, no significant plate transformation could be 
computed in order to further improve Chandra a strometry. 
Use of the USNO-B1 catalog ( IMonet et al.|2003[ ) yields one 
further possible coincidence (with a larger offset ~ l'-'l), 
which is of no help. 

On the optical side, we then re-computed the astromet- 
ric calibration of the large-field of view ACS image against 
the positions of 35 well-suited reference stars (i.e., not too 
faint, not saturated, not close to the CCD edges or to 
diffraction spikes) selected from 2MASS n identified in the 
mosaic of the two ACS chips. We measured the pixel co- 
ordinates of the 2MASS sources through gaussian fitting 
with the Graphical Astronomy and Image Analysis (GAIA) 
toorland we computed the pixel-to-sky coordinates trans- 
formation using the code ASTROtQ This yielded an rms of 
oy ps O'.'ll in our astrometric fit, accounting for the rms of 
the fit in the right ascension and declination components. 
Thanks to the pixel scale of the WFC, we neglected the un- 
certainty on the reference star centroids. Following |Lattanzi| 
et al. (1997), we also estimated the uncertainty in the reg- 
istration of the ACS image on the 2MASS reference frame. 
This is given as Otr — V3xo- re f/vN re f, where V3 accounts 
for the free parameters in the astrometric fit, o~ re f <> 0'.'2 is 
the mean positional error of the 2MASS coordinates, N re f 
is the number of 2MASS stars used for the astrometric cal- 
ibration. In our case, we obtain a tr = / /06. We finally con- 
sidered the C//015 uncertainty on the link of 2MASS to the 
International Celestial Reference Frame (ICRF). Thus, by 
adding in quadrature the rms of the astrometric fit and all 
the above uncertainties, we obtained that the overall posi- 
tional accuracy of our ACS' astrometry is Sr = (y.'13 (lc). 
By finally adding in quadrature this value to the 0'.'6 error 
on the Chandra coordinates we obtained an overall uncer- 
tainty of ~ 0'.'62 on the registration of the IE 1207.4-5209 
position on the ACS image. 

Results are shown in Figure [T] Source Z lies about 
l'/l and l'/3 off the Chandra HRC and ACIS position of 
IE 1207.4—5209, respectively. Such offsets are larger than 
the expected 99% accuracy of Chandra astrometry, which is 
estimated to be of 0'.'8 close to the centre of the field of view. 



As already concluded by Mignani et al. (2007a) and Wang 



et al. 



(20071, the association of source Z to IE 1207.4-5209 



based on positional coincidence seems very unlikely. 



3. Photometry 

Results from our HST test allow us to exclude any asso- 
ciation of source Z to IE 1207.4—5209. Thus, we can use 
the deep ACS images to set upper limits to any undetected 
source at the position of IE 1207.4—5209. We focus on the 
two deep ACS observations obtained through the F555W 
and F814W filters. To estimate count rates, we used a cir- 
cular aperture of (f.'5 ra dius and then we per formed aper- 
ture correction following Sirianni et al. (2005). Count rate 
to magnitude conversions was performed using standard 
ACS' photometric calibration provided by the data process- 



Mignani et aL 1 2007a I used an early relea se of the Guide 
Star Catalog 2 ( |Lasker et al.||2008|), whil e |Wang et al] ( |2007 l 
used the USNO-A2 Catalog ( |Monet|1998p . 

4 star- www.dur.ac.uk/~pdraper/gaia/gaia.html 

5 http://star-www.rl.ac.uk/Software/software.htm 



Table 1. Summary of the VLT/ ISAAC observations of the 
IE 1207.4-5209 field, with the observing epochs, the filter, 
the total integration time, the average seeing, and airmass 
values. 



yyyy-mm-dd 


Filter 


T(s) 


Seeing (") 


Airmass 


2003-01-27 


J 


1200 


0.83 


1.15 


2003-01-28 


H 


1000 


0.93 


1.14 


2003-02-11 


J 


2280 


0.93 


1.18 


2003-02-13 


Ks 


1000 


0.83 


1.32 


2003-02-16 


Ks 


3500 


1.17 


1.25 


2003-02-16 


H 


1750 


0.73 


1.30 


2003-02-19 


H 


1750 


0.77 


1.48 


2003-03-18 


H 


1750 


0.50 


1.44 



ing pipeline. Taking into account the observed background 
noise in a portion of the image surrounding the position of 
IE 1207.4—5209, we set a 3tr upper limit of nip 81 4w=28.1 
and niF555w=28.1. For completeness, we also computed the 
flux of source Z, which resulted to be mp S i4\y — 24.71±0.01 
and m i ?'555 W / = 26.80±0.05. Using the more recent WFPC2 
image, we computed m^ 814W = 24.83 ± 0.08, consistent 
with the ACS flux in the same band. The quoted values are 
not dereddened. 

Moreover, we include in our photometric study a se- 
ries of IR observations collected with the VLT and avail- 
able in the pu b lic ES Q arc hive (Programme 70. D-0436A). 
Pavlov et aL] ( |2004[ ) and |Fesen et all ( |2006[ ) reported 
preliminary results from such data. IR observations of 
IE 1207.4—5209 were performed in service mode with the 
Infrared Spectrometer And Array Camera (ISAAC) instru- 
ment at the VLT (Paranal Observatory) between January 
27th and March 18th 2003. The ISAAC Short Wavelenght 
(SW) camera was used, with a projected pixel size of 0'.'148 
and a field of view of 152x152 arcsec. Observations were 
performed through the J (A = 1.25//; AA = 0.29//), H(X = 
1.65//; AA = 0.30//) and K s (A = 2.16//; AA = 0.27//) 
band filters. To allow for the subtraction of the variable 
IR sky background, observations in each filter were split in 
sequences of shorther dithered exposures with integration 
times of 50 s in the H and K s bands and of 120 s in the J 
band along each point of the dithering pattern. The journal 
of observations is reported in Table 1. 

The total integration times over all nights were 3480 s 
(J), 6250 s (H), and 4500s (K s ). For each band, observa- 
tion were taken under photometric conditions with a seeing 
often better than l'.'O and airmass below 1.5. Atmospheric 
conditions were average with only the nights of February 
11th, 16th, and 19th affected by a humidity up to 40%. 
Twilight flat fields, dark frames, as well as images of stan- 
dard stars from the Persson et al. (1998) fields, were taken 
daily as part of the ISAAC calibration plan. We downloaded 
the data from the ESO public Science Data Archivqj and 
we reduced/calibrated them using the updated version of 
the ESO ISAAC pipeline^ For each exposure sequence, sin- 
gle frames were registered and coadded to produce a back- 
ground subtracted and cosmic-ray free image. We used the 
ACS 814W image as a relative reference frame to register 
the IE 1207.4-5209 position on the ISAAC images. 

e www.eso.org/sci/archive 

7 http://www.eso.org/observing/dfo/quality/ISAAC/pipeline 
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Table 2. Results of photometry close to the position of 
IE 1207.4-5209, based on HST and VLT data. Values are 
not corrected for interstellar reddening. 



Instrument/Filter 


Upper limit 


Source Z 


ACS /F555W 


28.0 


26.80 ± 0.05 


ACS /F814W 


28.1 


24.71 ±0.01 


ISAAC /J 


23.9 


21.53 ±0.07 


ISAAC /H 


22.7 


20.63 ±0.11 


ISAAC /Ks 


21.7 


20.53 ±0.18 



Photometry was performed using the SExtractor soft- 
wara^v2.4, which implements the "first moment" algorithm 
( Kron|1980 ). In view of the non-optimal sky conditions, we 
performed a photometric calibration of each image using a 
set of 25 2MASS stars as a reference. Our solutions turned 
out to be very good, with a r.m.s. of ~ 0.1 mag in the J, H 
and Ks bands, respectively. The resulting fluxes of source Z 
are J=21.53±0.07, H=20.63±0.11, Ks=20.53±0.18. Such 
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Fig. 2. Infrared-to- X-ray spectral energy distribution for 
IE 1207.4-5209. Flux limits from HST and VLT data -cor- 
rected for interstellar reddening assuming Nff = 1.3 x 10 21 



resul ts are in broad agre ement with Wang et al. 

with Fesen et al. (2006). The upper limits to trie emission 

of IE 1207.4-5209 areJ~ 23.9, H~ 22.7 and Ks- 21.7. 



4. Discussion 

Our test with HST, based on both absolute and relative as- 
trometry on multi-epoch images, firmly rules out any phys- 
ical association of Source Z with IE 1207.4—5209. 

What is source Z? Multicolor photometry, based on the 
HST and VLT datasets, points to an unrelated background 
red dwarf. Flux and colors are consistent with (although 
slightly redder than) an M5 star located at — 5 kpc, red- 
dened by E(B-V)~0.1. We note that the VLT upper limits 
based on our 2002 observations in the R band (R>27.1, 

with 



(I2OO7I) and cm a nd using the relations by Predehl & Schmitt ( 1995 ) 



and by Cardelli et al 



De Luca et al. 2004) are not consistent with the source 
flux as measured with HST in 2005. This could be due to 
some intrinsic variability of the dwarf star, and/or to con- 
fusion effects in the VLT images, due to the PSF wings of 
the two much brighter stars lying a few arcsec away to the 
South- West. Of course, in both cases the conclusion about 
the non-association of source Z with IE 1207.4—5209 would 
not change. 

Thus, IE 1207.4—5209 remains, as yet, unidentifi ed in 
the optical/IR, as all the other CCOs observed so far ( |Fesen 
et al.|[2006l |De Luca et~ld~1 120"0"8{ |Mignani et~aT]|2008 



20093 



3) 



111 

en 
.8, 



with the only possible exceptio n of the source 
in the" Vela Jr. SNR ( |Mignani et al.|2007b[ ). The upper lim- 
its to the optical/IR flux presented here are the deepest 
available so far for a member of the CCO class and corre- 
spond to an (unabsorbed) optical - to - X-ray flux ratio 
F8i4w/Fo.3~3keV ~ 5 x 10~ 6 . The spectral energy distri- 
bution for IE 1207.4-5209 is shown in Figure [2] 

Such upper limits virtually rule out the possibility of 
any stellar companion tied to the NS in a binary systerrF] 
At the distance of IE 1207.4—5209, taking reddening into 
account, only a very low-mass star [M — 0.1 M ) would be 
allowed. It is very unlikely that a binary system featuring 
such a low-mass star could survive the supernova explosion. 



(119891) -are s hown, togeth e r wit h 
the limits from Spitzer data given by Wang et al. (2007). 
The expected flux from a fallback disk, assuming i = 60 , 
B= 8 x 10 10 G, m = 8 x 10 11 g s~\ r out = 2 x 10 12 cm 
(see text), consistent with the upper limits, is plotted as a 
red dash ed line. The X-ray sp ectrum as observed by XMM- 
Newton ( De Luca et al.|2004[ ) is also plotted as a solid black 
line. 



Combining the constraints from the optical/IR flux lim- 
its to t he ones derived from the t iming properties of the 
source ( [Gotthelf fc Halpern 2007), we may test the hy- 
pothesis ^oFTh^ - exisTelice - ot aclisk formed of supernova 
fallback material, surrounding the NS. Such a debris disk 
has been invoked as a possible explanation of several puz- 
zling properties of the X-ray source. In order to compute 
the expected optical/IR flux from such a putative disk, 
we consider a geometrically thin, optically thick fallback 
disk, locally emitting as a blackbody. Follo wing a stan- 
dard approach, we included two contributions (Perna et al. 
2000): (i) viscous dissipation i n the disk, yielding a teir T 
perature profile T(r) oc r~ 3 / 4 (Shakura & Sunyaev 1973); 
(ii) reprocessing of X-rays emitted by the central object. 
Assuming Lx = 2.2 x 10 33 erg s _1 for a distance of 2.2 kpc 
( De Luca et al.|[2004] ) and a standard disk structure yield s 
T(r) = 1100if(l - ?7) 2 / 7 (i? Q /r) 3 / 7 ( [Vrtilek et alT||l990| ), 
where 77 is the X-ray albedo of the disk. The total flux is 
obtained by integrating the emissivity over the disk surface 
and taking into account the inclination of the disk with re- 
spect to the line of sight, as well as the source distance. The 
inner radius of the disk is set at the distance from the star 
where the magnetic pressure of the rotating dipole of the 
neutron star disrupts the disk itself. Taking into account the 
role of viscosity in the disk, such a radius is tm ~ 0-5Rm 
where Rm is the so-called Alfven radius (e.g. Frank et al. 



http://www.astromatic.net/ 
9 The limits on P leave some room for a long period system 
featuring a small star. 



2002). The outer radius of the disk may be constrained by 
the flux limits. We assumed a disk inclination i = 60°. The 
contribution of the reprocessed X-ray flux critically depends 
on the poorly known value for the X-ray albedo 7/. Although 
a val ue r\ = 0.5 has been assumed in severa l investigations 
(e.g. Vrtilek et al. 1990 Perna et al. 2000), a much larger 
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value r\ — 0.97 was evaluated by Wang et al. (12006), based 
on the detection of a disk around AXP4U 0142+61. Thus, 
we used r\ = 0.97 (yielding a much lower disk luminosity) 
in our study. 

In the above assumptions, we computed the expected 
flux in a specific filter band as a function of the NS mag- 
netic field and disk accretion rate m, for a set of values of 
the disk outer radius. We repeated such exercise for all of 
the bands of our optical/infrared dataset. For completeness, 
we also used the Spit zer 4.5/xm and 8.4/im bands studied 
by Wang et al. (2007). The most constraining limit turned 
out to be those lrilhe HST F555W and F814W bands re- 
ported here. In the allowed range of parameters, the model 
is rather insensitive to the value of the disk outer radius. 
The flux contributions are almost negligible in all bands 
for r out larger than ~3x 10 12 cm. Results are plotted in 
Figure [3] where a red line marks the rh-B field region al- 
lowed by the combination of the F555W and F814W flux 
limits. The region ruled out by HST limits is colored in yel- 
low. For comparison, we also marked (in orange) the region 
formerly ruled out by the less constraining Spitzcr/IRAC 
4.5/im limit. Flux limits in other bands are less constraining 
than the Spitzer one. 

In order to use the tim ing constraints (P < 2.5 x 10~ 16 
Gotthelf & Halpern 2007 ) , we note that the interaction of 



the putative disk with the rotating neutron star magneto- 
sphere should yield different regimes of angular momentum 
transfer, according to the relative positions of Tm with re- 
spect to the neutr on star light cylinder as we ll as to the 
corotation radius (Illarionov & Sunyaev 1975). We evalu- 



ated the expected neutron star P as a function of the NS 
magnetic field and disk rh in such different regimes (the so- 
called ejector, propeller and accretor regimes) using stan- 



dard relations (e.g. |Menou et al.||1999| |Zavlin et al.|2004 ). 
This allows to identify a rh-B field region allowed by the 
existing limit on P, which is overplotted in Figure wl The 
line marking the P limit has been drawn arbitrarily (in 
dashed style) in the small region connecting the propeller 
regime to the ejector regime (where standard relations for 
the propeller torque are not valid) , to visualize the reduced 
efficiency of the propeller effect as the magnetospheric ra- 
dius approaches the light cylinder. The region ruled out by 
X-ray timing is colored in grey. A region in the accretor 
regime yielding a luminosity larger than 2 x 10 33 erg s _1 
(i.e. the total X-ray luminosity of IE 1207.4—5209 assum- 
ing a distance of 2 kpc) is also ruled out. It is marked in 
red. 

The region allowed by both the P limit and the opti- 
cal/infrared flux limits is coloured in green in Figure [3] We 
can infer some interesting indications about the possible 
role of a fallback disk. In the absence of any model describ- 
ing a quantitative relation between the accretion rate and 
the phenomenology of the spectral features, we will focus 
here on the issue of the luminosity of the hot spot. The 
maximum allowed disk rh is slightly less than 10 12 gs' 1 in 
the propeller regime. However, as a consequence of the cen- 
trifugal barrier, only a minor fraction (if any) of such mass 
rate would enter the Alfven radius. On the other hand, 
in the accretor regime, which is allowed only for very low 
magnetic field values (smaller than a few 10 8 G !), we would 
expect an accretion rate smaller than ~5x 10 10 g s _1 . In 
both the propeller and the accretor pictures, the accretion 
rate on the INS could yield a luminosity of ~ 10 31 erg s _1 . 
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Fig. 3. Limits on the neutron star magnetic field and rh 
of a put ative fallback disk in IE 12 07.4-5209, set by X-ray 
timing ( |Gotthelf fc Ha lpern 2007) and by optical/infrared 
flux limits. Golors are as follows: (grey) - ruled out by X-ray 
timing; (red) - ruled out by total X-ray luminosity (assum- 
ing d=2 kpc); (orange) - consistent with X-ra y constraints 
but a lready ruled out by former Spitzer limits (Wang et al. 
1200 7|; (yellow) - consistent with X-ray constraints but ruled 
out by our HST limits; (green) - parameter region consis- 
tent with all constraints. 



This is about two order of magnitude smaller than the X- 
ray luminosity of the hot thermal component (~9x 1 32 
erg s -1 ) seen in t he X-ray spectrum of IE 1207.4—5209 (De 
Luca et al.(2004| ) 



Of course, such results are model-dependent and should 
be taken with some caution. However, they point to an al- 
ternative explanation for the origin of the high temperature 
and luminosity of the hot spot on IE 1207.4—5209. A very 
similar picture emerged for PSR J1852+0040, the CCO at 
the center of t he Kes 79 supernova remnant. |Halpern fc| 
Gotthelf (2010) succeded in measuring the P of the INS 
and ruled out fallback accretion as a heating mechanism 
for the "hot spot" seen in that source (such hot regions 
are ubiquitous among CCOs). The same authors concluded 
that other possible explanations of the thermal anisotropics 
in CCOs (e.g. localized magnetic heating in a sunspot-likc 
configuration, or anisotropic conduction through a strong 
poloidal magnetic field in the INS) also suffer of several dif- 
ficulties, which also apply to the case of IE 1207.4—5209. 

It would also be interesting to constrain the total mass 
contained in the putative fallback d isk, in view of the new 
flux limits. However, as discussed by Beckwith et al.|( |1990[ ), 
the disk is expected to be optically thick at wavelengths 
shorter than ~ 1 mm. Thus, a measure of the disk emission 
at very low frequenc ies would be requ ired to get a direct 
mass estimate. While Wang et al. ( 2006 ) , in the case of AXP 
4U 0142+61, could extrapolate the observed spectral shape 
of the disk to set a flux limit in the mm band, our upper 
limits do not allow us to set constraints to the emission of 
the system at very low frequencies, hampering any direct 
mass estimate. On the other hand, we can use the limits on 
the disc m, coupled to a model for fallback disk evolution 
as a f unction of tim e (Cannizz o et al. 1990 Menou et al 



1999 Chatterjee et al. 



2000{ l De Luca et al.||2006p . Such 
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models expect that, after an initial (uncertain) phase during 
which supernova fallback material settles into a disk, the 
disk evolution is ruled by viscous processes and obeys a 
law rn ex t _Q , with a in the range 1.18-1.25 (depending 
on the disk opacity) . Of course, the disc interacts with the 
rotating, magnetized INS and thus affects its spin-down 
history. We assumed an age of 10 4 yr for the system, a 
current disc m = 5x 10 11 g s _1 , a dipole B=8 x 10 10 G field 
for the neutron star (consistent with a propeller regime, see 
FigJ3J), and a dec ay index a = 1.25. Th e model (similar to 
the one used by |De Luca et al. 2006 ) allows to estimate 
the maximum disk mass at the onset of the disk viscous 
evolution (t = to, assuming the INS period P(to) = 424 
ms) consistent with all of the flux and timing limits. Such 
an "initial" mass turns out to be as small as 10 -6 M Q . The 
current disk mass (at t = t + 10 4 yr), according to the 
model, would be about 10 4 times smaller. 



5. Conclusion 

Relative and absolute astrometry on HST datasets rules out 
the association of a faint optical source to IE 1207.4—5209. 
Using HST and VLT data, we derive very deep optical and 
infrared upper limits to any emission from IE 1207.4—5209. 
Our optical/infrared upper limits, coupled to the con- 
straints derived from the limits on P, argue against the 
possibility that accretion from a disk of fallback material 
could explain the existence of a luminous, hot thermal com- 
ponent in the X-ray spectrum of this source. As for other 
CCOs, the large surface thermal anisotropy remains unex- 
plained. Using published models for the evolution of fall- 
back disks as a function of time, we also estimated that the 
mass of any disk currently surrounding the neutron star (at 
an age of ~ 10 4 yr) should be smaller than ~ 10 -10 M Q . 

The most likely picture for IE 1207.4-5209 is that of a 
young, isolated neutron star, born with a rotation period 
very similar to the current one, spinning down at a tiny rate, 
ruled by a dipole field smaller than ~ 3x 10 11 G. Most likely, 
the peculiar cyclotron features seen in the X-ray spectrum 
of the source are generated very close to the star surface. 
A very low level of accretion (< 5 x 10 10 g s~ 4 ) cannot be 
excluded, and its role for the production of such features 
should be explored. Using the current generation of ground- 
based or space observatories, more sensitive optical/IR data 
would require a very large investment of observing time. 
The James Webb Space Telescope (JWST), to be launched 
in 2014, will provide the required sensitivity and spatial 
resolution in the near-IR to improve on our limits. 
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